Electron linewidths of wide-gap insulators: excitonic effects in LiF 
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Based on a recent exchange-correlation kernel developed within Time-Dependent-Density- 
Functional Theory we derive a practical and general expression for the three-point vertex function. 
We show that excitonic effects in LiF strongly modifies the low-energy electron linewidths leading 
to linear scaling with quasiparticle energy. We also prove that, in contrast to previous results for 
the electron gas, simple metals and semiconductors, vertex corrections in the self-energy and in the 
screening function do not compensate each other. 
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The experimental quasiparticle band structure of bulk 
metal and semiconductor systems has been successfully 
explained by the GW self-energy scheme Q, 0, in its 
simplest non self-consistent GqWq implementation. Sim- 
ilarly, the quasiparticle linewidths of simple and noble 
metals have been studied extensively 0, but a first- 
principles description of the electronic contribution to 
the electron/hole linewidths in semiconductor and insu- 
lators is not yet available. The reason is that the low- 
energy quasiparticle dynamics in semiconductors tends 
to be dominated by inelastic phonon scattering, the elec- 
tronic contribution playing a minor role. However this 
scenario changes drastically when the quasiparticle en- 
ergy is larger than the minimum energy required to ex- 
cite an electron-hole pair. Above this threshold (that is 
zero in metals and approximatively twice the band gap 
in insulators) the rapidly increase in density of electron- 
hole pairs dominates the quasiparticle damping, resulting 
in a mainly electronic contribution to the lifetime. It is 
well known that in insulators, at difference with metals, 
the attractive interaction between electrons and holes can 
lead to the formation of a bosonic-like excitonic state Q . 
Excitons modify remarkably the optical and energy-loss 
spectra and, consequently, the microscopical mechanisms 
responsible for the quasiparticle damping. This effect is 
stronger in wide-gap insulators like LiF. 

In this communication we tackle the problem of evalu- 
ating the impact of the excitonic effects on the quasi- 
particle (QP) dynamics of LiF, using a simplified ver- 
tex function in the electronic self-energy. An efficient 
approximation for the three-point many-body vertex 
function is given in terms of the two-point exchange- 
correlation kernel f xc recently developed in the 
framework of Time-Dependent-Density-Functional The- 
ory (TDDFT) . As result we show that the electronic 
linewidths of LiF display a linear dependence as function 
of the QP energy, that can be traced back to the incip- 
ient excitonic effects induced by f xc 0. Understanding 
the inelastic mechanism which dominates the phase co- 
herence time is crucial to the field of quantum transport 
in mesoscopic and nanostructred materials. Thus this 



work is the first steps toward a full first-principles de- 
scription of the quasiparticle dynamics of semiconductor 
and insulators. 

In the usual one-shot GqWq self-energy scheme, it is 
assumed that a basic Density-Functional-Theory (DFT) 
calculation ?] provides good approximation for QP wave- 
functions and electronic screening (dominated by col- 
lective excitations, plasmons build from independent 
electron-hole transitions, i.e. excitonic effects in the 
screened Coulomb potential Wq are neglected). The 
QP lifetime r, can be calculated with the Fermi golden- 
rule, using this non interacting Wo as scattering poten- 
tial: r" 1 = -2 \Q if | 2 Im [W (E, - E f )], where \i), 
|/) are the initial and final state, with Wq matrix ele- 
ments Qj/, and energies Ei Ef such that E{ — Ef > 0. 
This scattering scheme, also known as "on-mass shell" 
approximation to the GqWq linewidths, provided valu- 
able insight into the electron/hole linewidths of met- 
als 0. Therefore, in Fig. we estimate the electronic 
line widths of LiF (boxes) within this approximation. As 
Wo is calculated in terms of non-interacting electron- 
hole pairs rf 1 — when E t — Ef < E gap (E gap the DFT 
gap): quasiparticle states with energy Ei < 2E gap have 
zero line width (infinite lifetime). These states are indi- 
cated by the dashed area in Fig. Above this region a 
quadratic energy dependence of the line width is recov- 
ered, as in metals g|. 

As the short-range screened Coulomb repulsion mod- 
ifies drastically the polarization function in LiF, one is 
tempted to apply the previous "on-mass shell" scheme 
to analyze the role of excitonic effects on the quasipar- 
ticle dynamics. This would correspond to replace Wo 
by the screened coulomb potential W obtained from the 
many- body Bethe-Salpeter equation (BSE) [1HE3]- In 
practice, the BSE sums all the possible binary collisions 
between electrons and holes, providing a consistent and 
successful framework for the calculation of the interacting 
polarization function. However, the BSE is computation- 
ally very demanding and it becomes unpractical when the 
microscopical dielectric matrix e(q, w) must be calculated 
for a large set of transfer momenta q and frequencies uj, as 
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it is the case for the calculation of linewidths Q . To by- 
pass this difficulty we compute efq, us) within a TDDFT 
framework, using an f xc kernel |5| that mimics well the 
BSE results E3- This performance is illustrated in Fig. [21 
for the loss function e _1 (q, lj), that is the relevant quan- 
tity to build the screened Coulomb potential W. From 
Ref. we know that TDDFT reproduces the experimen- 
tal loss function, therefore comparing the TDDFT and 
random phase approximation (RPA) results of Fig. we 
see that RPA misses the strong weight of the loss spectra 
just above the band-gap. Consequently the inclusion of 
excitonic effects in this GqW calculation translates into a 
drastic change of the quasiparticle decaying rates (red cir- 
cles in Fig. ^| compared to the RPA results (blue boxes) . 
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FIG. 1: (Color on-line) Left panel: calculated DFT band- 
structure of LiF (here Ecbm and E ga p stand for the DFT con- 
duction band minimum energy and the energy gap). Right panel: 
Electron linewidths calculated "on mass— shell" as function of the 
single-particle energy. Boxes: RPA Go Wo- Circles: TDDFT- 
based vertex correction to the self-energy, i.e, a GqWT^P ddft 
approach that turn out to be very close to a simpler GqW calcu- 
lation (see text). The dashed area denotes the forbidden energy 
region for quasiparticle decay into electron-hole pairs. Error bars 
represent the the oretical uncertain due to the "zero-broadening" 
extrapolation llll . 



This simple scattering approach, though appealing, 
lacks of theoretical consistency. Exchange-correlation ef- 
fects have been included only in the polarization function, 
while, in the spirit of the original work of Hedin 0, they 
should be included in the self-energy as well. However we 
will show below that the results obtained within a proper 
treatment of self-energy and polarization effects do not 
deviate appreciably from the previous GqW results. We 
start the derivation from the definition of the self-energy 
operator £ (1, 2), given by [t| 

E(l,2)=z J d34W(l+3)G(l,4)r(4,2;3). (1) 

Here G(l,2) is the interacting Green's function and 
r (1, 2; 3) the irreducible vertex function (numbers stands 
for space, time and spin coordinates). The screened 



Coulomb interaction W is: W (1,2) = v(l,2) + 
J d34u(l,3)x(3,4)W(4,2) where i?(l,2) is the bare 
Coulomb interaction and \ the irreducible polarization 
function: 

X(l,2) = -i[ d34G(l,3)G(4,l)f (3, 4; 2). (2) 

Thus, given an approximation for T the self-energy is 
completely defined trough Eqs. Q|2J) plus the Dyson 
equation for G. Electron-hole effects are embodied in the 
vertex function T that appears in the self-energy directly, 
in Eq. QJ, and trough the polarization function, Eq. J3J). 
The interplay between those two effects has been strongly 
debated in the last years, using different approximations 
for r, and different levels of self-consistency in the solu- 
tion of Dyson equation. However all the systems analyzed 
in the past are characterized by moderate, if not absent, 
excitonic effects in the polarization function. Thus even 
if the use of two-pointy DFT-based [HE! El or finite 
order vertex functions F can be justified in the case of 
the homogeneous electron gas or simple semiconductors, 
they will be inadequate in the case of wide-gap insula- 
tors (e.g. LiF), as well as in the case of other strongly 
correlated systems. Next we derive a TDDFT approx- 
imation to the vertex function T following the spirit of 
Ref. [5j : to reproduce the diagrammatic expansion of T 
obtained within Many-body perturbation theory. 
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FIG. 2: Calculated loss function of LiF for momentum transfer 
q along the direction Y U . Continuous line: TDDFT calculation 
with an f xc kernel that mimics excitonic effects [jj. Dashed line: 
RPA. 

In the non self-consistent scheme the BSE expresses T 
in terms of the independent particle Wq (calculated using 
Eq. © assuming f (1, 2; 3) = 6 (1, 2) 5 (1, 3)) and the bare 
DFT Green's function Go as: 

f (1,2; 3) = 5 (1,2) S (1,3) + 

iW (1, 2) / d67 G (1, 6) G (7, 2) f (6, 7; 3) . (3) 
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When this vertex T is inserted in Eq. J3J), the correspond- 
ing equation for \ correctly describes excitonic effects in 
the polarization function at the BSE level |3| • Even if T 
is an highly non-local, three-point function, it has been 
recently shown that, as long as we are interested in the 
two-point polarization function %, Eq. can be cast in 
terms of the two-point exchange-correlation kernel f xc 
of TDDFT HQ: 

X (1, 2) = xo (1, 2) + / d34 Xo (1, 3) f xc (3, 4) X (4, 2) . 



(4) 

Here xo (1, 2) = -iG (1, 2) G (2, 1) gives the DFT po- 
larization function. At this point if we take the exchange- 
correlation potential corresponding to f xc as a local ap- 
proximation to the self-energy, then the vertex func- 
tion can be easily contracted into a two-point function: 
f (6, 7; 3) = f loc (6, 3) 6 (6, 7) M , with 



(1,2) 



~] -l 



6(1,2)- / d3/ xc (l,3) X o(3,2) 



(5) 



the 



Thus Eq.JTJ gives 2(1,2) 

iW TDDFT(i+^ 2 ) G (1,2), in terms of 
TDDFT effective potential W TDDFT (1, 2) 
/ d3 v(l, 3) [S (3, 2) - / d4 (v (3, 4) + f xc (3, 4)) X o (4, 2)] "\ 
From this self-energy the lifetime of a generic conduction 
state c with momentum k is given by 

7--^ = -2ft- 1 J2 E^'( k q G i)^'( k q G 2) 

G!,G 2 q,c' 

Im[W^ T {^-e c ,^)}, (6) 

with p„„/ (kqG) = (nk|e 4 ' q+G ^' r |n'k — q), G a reciprocal 
space vector and VL the crystal volume. Different expres- 
sion for Fioc, based either on local-field factor of the ho- 
mogeneous electron gas |l3j or on time-dependent local- 
density approximation (TDLDA) 0, Q] have shown 
that inclusion of local vertex corrections in both S and x 
almost cancel out, i.e., T[ oc in E undresses the exchange- 
correlation effects included in the polarization function 
X- However such approximations for f xc produce optical 
spectra very similar to RPA, in disagreement with exper- 
iments. This important drawback of a TDLDA f xc has 
been recently related to the long-range nature of the ker- 
nel, f xc (r, r'; ui) ~ — a (u>) /|r — r'| that partially coun- 
teracts the repulsive Hartree contribution 0, . The 
stronger the electron-hole effects are, the larger is the 
correction embodied in a. In the case of wide-gap in- 
sulators like LiF there is a large region of frequencies 
and transfer momenta q where f xc is stronger than the 
Hartree term (i.e. a > 1). This leads to unphysical 
linewidths: for a large energy range Im (2), and hence 
t -1 , has a wrong sign ! This result is visualized by notic- 
ing that with respect to a GqWq calculation a change of 



sign of r 1 is controlled by the sign (v + f xc ) , that is pro- 
portional to (1 — ct(u>)). A similar result was obtained 
in Ref. [HI looking at the high q limit of the TDLDA 
kernel that goes as f xc ~ q 2 . The reason for this impor- 
tant failure of a two-point vertex function is connected 
to the imposed reduction of the non-locality from the 
original, three-point vertex function. In physical terms 
Tioc overestimates the intensity of the vertex correction 
because two incoming particles (entering in 1 and 2 in the 
exact vertex function T) are supposed to coexist at the 
same time-space point. To overcome this difficulty we 
decided to release the constrain on the spatial locality 
and define a TDDFT vertex function Ttddft (1, 2; 3) 
such that, for a given / xc (l,2), Ttddft is consistent 
with Eqs. To this end we recall that in Ref. |5| we 

derived a diagrammatic expression for f xc in terms of the 
screened coulomb potential Wo, that to first order reads: 
fxc = Xo 1 X Xo-L with x' 1 ' the first order expansion 
of Eq. (PJ) in Wo |L7j] • From this f xc we get an approxi- 
mation for the vertex function T^ DDFT (1, 2; 3) imposing 
that once plugged in Eq. (01 it reproduces Eq. 10} for x- 
By inspecting Eqs. we obtain 



Vddft (1»2; 3) = S (1,2) 8 (2, 3) + 

iWo (1, 2) / d4 Go (1, 4) Go (4, 2) f loc (4, 3) , (7) 



It is crucial to observe that T^ DDFT is not a first or- 
der vertex, as r; oc sums an infinite number of diagrams. 
Eq. J7J) can be easily generalized to give higher order ap- 
proximations for r, consistent with the high order cor- 
rections to f xc of Ref. 5]. As it is commonly done we 
neglect dynamical effects in the BSE , i- e - we assume 
W (1,2) w W (ri,r 2 ;w = 0) in Eq. 0. This approxi- 
mation is motivated in the present case, as we are in- 
terested in the low-energy electronic linewidths neglect- 
ing self-consistency effects. We have verified numerically 
that for LiF, Si, diamond and SiC>2 Ttddft 0-> 2; 3) is an 
excellent approximation to the "true" BSE vertex func- 
tion r. 

Now we can study the quasielectron lifetime in this 
approximation for the vertex function and for the elec- 
tronic self-energy. To do so, we use as above the "on 
mass-shell" approximation, i.e., the lifetime is given by 
the imaginary part of 2 = GoI4 7 r T 1 ^ )I)FT evaluated at 
the DFT energies: 

T ck = - Im [( ck |2 (r, r'; e ck ) + E (r, r'; e ck ) |ck)] . (8) 

Here the linewidths are computed as an average of the 
"left" and "right" self-energies, 2 and 2 in order to re- 
store the proper r, r' symmetry of the self-energy. Using 
Eq. J7J r^ 1 can be simplified by performing the energy 
integration in the complex plane and exploiting the pole 
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structure of T 



TT 1 = T.- 1 



(1) 

TDDFT' 



ck — 'ck,0 

-20- 1 x: E^[^f G T T (q.^-e c . k - q )] 

Gi,G 2 q,C 

i?e [(T£ (kqGi) + O (kqd)) (kqG 2 )] , (9) 

where t~£~ corresponds to the standard GqW approxi- 
mation, and 

GiG 2 Q c 2 v 2 

Pdv 2 ( k - qQG 2 ) Pc 2 v 2 (k - QqG) 

(e c k — e c 'k- q — e C2 k-Q + e« 2 k-q-Q) 1 , (10) 



an interpretation of the quasiparticle scattering based 
only on independent-particle processes losses meaning. 
The results of the present work allow for the systematic 
analysis of the role of excitons in quasiparticle excita- 
tions and response functions of extended and low dimen- 
sional systems, where the standard GqWq approximation 
fails HQ. 
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r-,(kqG) = «- 1 E^^iH^^g.g 

G1G2Q C 2 -"2 

p* c , C2 (k-qQG 2 )p„ 2C2 (k-QqG) 

(e c k — e c 'k-q — e« 2 k-Q + e C2 k-q-Q) 1 ■ (11) 

Eq. constitutes the main basic result of this commu- 
nication, and can be easily extended to the quasihole 
linewidths. Eq. (jSJ must be compared with Eq. (0- In 
the case of weakly interacting systems the two equations 
with a TDLDA f xc give very similar quasiparticle correc- 
tions to the gap and electron linewidths |l2lll3lll4| . But, 
as short-range correlations become important Eq. © 
tends to give non-sensible results (negative linewidths) 
because of the wrong sign of W TDDFT . In Eq. (0, in- 
stead, the term (T™, + T^, ) reflects the spatial non- 
locality of ^j-ddft strongly reducing the weight of 
W TDDFT . Consequently the final expression for r"^ 1 is 
given by t^ q plus a small vertex correction that does 
not change appreciably the results of a simpler GqW cal- 
culation. The fundamental practical result of this work 
corresponds to the solution of Eq. 10 for LiF, shown in 
Fig- HI The overall effect of excitons in the linewidths is 
huge [l^l : the linewidths up to 3 eV above the forbidden 
region display a linear dependence with energy while the 
RPA are almost zero because of the slow rise of the RPA 
loss function (see Fig[5J). A similar energy dependence 
has been observed in highly correlated materials ||- In- 
stead the present linear dependence of the linewidths is 
due to the combination of an almost constant density- 
of-states close to the conduction band minimum and 
to a "step-like" energy dependence of loss function (see 
FigJ2J. Furthermore, the quasiparticle linewidths are not 
exactly zero in a small energy window of 0.5 eV in the 
forbidden region. This effect can be traced back to the 
excitonic-induced transfer of oscillator strength in the 
dynamical dielectric function below the gap. This result 
is consistent with the fact that exciton dynamics is dic- 
tated by vertex-correction to the self-energy, therefore 
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